Can mesoscopic fluctuations reverse the supercurrent through a disordered Josephson 

junction? 
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We calculate the Josephson coupling energy U.j{(f)) (related to the supercurrent / = {2e/K)dU.j /d<j)) 
for a disordered normal metal between two superconductors with a phase difference 4>. We demon- 
strate that the symmetry of the scattering matrix of non-interacting quasiparticles in zero magnetic 
field implies that Uj{(j>) has a minimum at = 0. A maximum (that would lead to a yr-junction or 
negative superfluid density) is excluded for any realization of the disorder. 

PACS numbers: 74.80.Fp, 72.f5.Rn, 73.63.Rt, 74.40. -Hk 



The question in thje title was posed ten years ago by 
Spivak and KivelsonJJ in their search for mechanisms that 
would lead to a negative local superfluid density in a dirty 
superconductor. A negative instead of a positive super- 
fluid density means that the Josephson coupling energy 
Uj{(f)) has a maximum instead of a minimum for zero 
phase difference (f> of the superconducting order param- 
eter. The supercurrent I{4>) = {2e/h)dUj /dip then has 
the opposite sign as usual for small 4>. This has a vari- 
ety of observable consequences, including a ground state 
with a non-zero supercurrent J3 Aharonov-Bohm oscillcy-i 
tions with period h/Ae, and negative magnetoresistance.El 
A Josephson junction with a negative superfluid density 
is known as a 7r-junctionp because the ground state for 
large magnetic inductance is close to (j) ~ n instead of as 
usual at (j) — 0. 

The kjMwn mechanisra for the creation of a tt^ 
junctiontl Q in equilibriumO in a non-magnetic materialcl 
requires strong Coulomb repulsion to create a local- 
ized spin on a resonant impurity level. Spivak and 
Kivelson asked the question whether purely one-electron 
conductance fluctuations might be sufficient to pro- 
duce a locally negative superfluid density near the 
insulating state. A suggestive argument that this 
might be possible comes from the relative magni- 
tude of the mesoscopic sample-to-sample fluctuations 
of the supercurrent in a disordered supercpnductor- 
normal-metal-superconductor (SNS) junction.Ll The ra- 
tio {SPy^"^ / {I) ~ e^/hG of the root-mean-squared fluc- 
tuations over the mean supercurrent is <C 1 if the con- 
ductance G of the normal metal is large compared to 
the conductance quantum /h. The ratio becomes of 
order unity on approaching the insulating state, suggest- 
ing that the supercurrent might have a negative value in 
some samples. 

Of course, the root-mean-square amplitude of the su- 
percurrent fluctuations does not distinguish between a 
positive and negative sign, so that this argument is only 
suggestive. We were motivated to settle this issue be- 
cause of recent experiments on localization in quasi-one- 
dimensional superconductors.! This has renewed the in- 
terest in the fundamental question whether mesoscopic 
fluctuations are sufficient or not to create a negative lo- 



cal superfluid density. The answer, as we will show, is 
that they are not. 

The two superconductors that form the SNS junction 
have order parameters Ae*"^/^ and Ae"**^/^. The contacts 
to the normal metal have N propagating modes at the 
Fermi energy Ep^ so that the elastic scattering by the 
normal metal at energy E = Ep -\-£ is characterized by a 
2N X 2N scattering matrix 5(e). The two properties of 
S that we use are that it is analytic in the upper half of 
the complex e plane and that it is a symmetric matrix, 
5(e) — 5(e)^, when time-reversal symmetry is preserved. 

The starting point of our calculation is the relationship 
derived in Ref. ^ between the Josephson coupling energy 
Uj{(j)) in equilibrium at temperature T and the scattering 
matrix, 

oo 

Uj = -2kBT J2 In Det [f - 5A(iw„)5Ar(ia;„)] . (f ) 

n=0 

The summation runs over the Matsubara frequencies 
ujn = (2n -I- l)nkBT. The AN x AN scattering matrix 
5jv(e) describes the elastic scattering by disorder in the 
normal metal of non-interacting electron and hole quasi- 
particles with excitation energy e, 



5Ar(e) 



(V 





5(-e) 



(2) 



The analytical continuation of 5jv from real to imagi- 
nary energy (e —>■ iu>) follows from 5(e) S{iuj) and 
5(— e)* — > S{iuj)*. Similarly, the matrix 5^1 (e) describes 
the Andreev reflection from the superconductors. 



Sa{£) = a{e) ,A0/2 







iA0/2 



a{e) = e/A - ly/l-eyA^. 
Here A is a 2N x 2N diagonal matrix with elements Aj 



(3a) 
(3b) 



f for 1 < j < iV and A,- 



-1 for iV-l- 1 < j < 2N. 



Eq. (yj) differs from the usual representation of the 
Josephson energy as a sum over the discrete spectrum 
(e < A) plus an integration over the continuous spec- 
trum (e > A). Each term in the sum (|l]) is combined 
from contributions out of the discrete and the continuous 



f 



spectrum. We will now show that each of these combi- 
nations is minimal for (j) = 0, although the contributions 
from the discrete and continuous spectrum separately are 
not. 

Let us abbreviate 



(4) 



Using the identity In Det = Tr In in Eq. (|l|) one can 
calculate the first and second derivative with respect to 
(j) of the energy Uj{(f>). The first derivative is given by 



dUj 



ImTr [hii{u!n) - /l22(Wn)] 



(5) 



where the N x N matrices hu and /i22 are blocks of the 
matrix 



H = Z*Z{1 + Z*Z)-' = ( J" ^2 
^ ' * n-21 n-22 



(6) 



The first derivative is equal to the supercurrent and van- 
ishes at (/) = as dictated by time-reversal symmetry. 
For the second derivative we obtain 



^AkeTY, ReTr [/i2K)/2*iK) 

ri=0 

+ /ll2(Wn)/l2l(Wn)] , 

fll fl2 
/2I /22 



F ^ Z{1 + Z*Z)-'^ 



(7) 

(8) 



At (f) — the symmetry of 5* implies that F = and 
H — H\ hence — fl2, /121 — /ii2- Therefore, every 
term in the sum ^) is positive. We conclude that the 
Josephson energy Uj{(j)) has a minimum at = 0: 



dUj 



0, 



d^Uj 



6=0 



> 0. 



(9) 



6=0 



This concludes the proof that mesoscopic fluctuations 
can not invert the stability of the SNS junction at zero 
phase, excluding a mechanism foi: the creation of a tt- 
junction proposed ten years agoE The proof holds for 
non-interacting quasiparticles in zero magnetic field at 
arbitrary temperature, for any disorder potential and any 
dimensionality of the junction. As a final remark, we 
conjecture (and have a proof for N — 1) that Eq. (1) 
implies dUj/d(j) > in the entire interval < < tt in 
the presence of time-reversal symmetry. 

We thank Boris Spivak for urging us to solve this prob- 
lem and for valuable discussions. This work was sup- 
ported by the Dutch Science Foundation NWO/FOM. 
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